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The cross section for massive fermion production in two-photon collisions was examined at next- 
to- leading order in QCD/QED for general photon helicity. The delta function (virtual+soft) part 
of the differential cross section was analytically integrated over the final state phase space. Series 
expansions for the complete differential and total cross sections were given up to tenth order in 
the parameter (3. These were shown to be of practical use and revealed much structure. Accurate 
parametrizations of the total cross sections were given, valid up to higher energies. The above results 
were applied to top quark production in the region not too far above threshold. The cross section 
was shown to be quite sensitive to a 3 in the appropriate energy region. 
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I. INTRODUCTION 

m 

High energy photons may be produced by backscattering laser light off high energy e~ or e + beams. In addition, 
high degrees of polarization are possible and the photons may carry a large fraction of the electron energy. Photon- 
£C) • photon collisions also arise naturally as a background in e + e~ collisions. One major motivation for constructing a 
77 interaction region at a high energy next linear collider (NLC) is to produce Higgs bosons on resonance via 77 
fusion, which also allows direct determination of the if 77 coupling, which is sensitive to possible non Standard Model 
charged particles of large mass that may enter in the triangle loop. Using polarized photons allows one to control 
the backgrounds arising from 77 — > bb for an intermediate mass Higgs [Q . This background has now been studied 
; including QCD @-| and electroweak l6f| corrections. 

In this paper we will consider in some detail the process 77 —fff + X in the region not too far above threshold, 
making use of the complete analytical results presented in [Q] , which include photon polarization. We will demonstrate 
the usefulness of 77 — > tt + X in determining a s precisely. We extend the analytical results presented in Q by 
1 integrating and obtaining analytical results for the single integral (virtual+soft) part and by series expanding the entire 
differential and integrated cross section to order /3 10 , where j3 is the massive fermion velocity in the soft radiation 
limit. Such an expansion is shown to be of practical use, not too far above threshold, and it demonstrates many 
interesting features of the corrected cross sections. We have also provided parametrizations of the total integrated 
cross sections valid up to higher energies. 

As the diagrams involve only QED-like vertices, the process under consideration is quite fundamental in nature. 
The fact that complete analytical results have been absent, until recently, reflects the lack of experimental feasibility of 
directly colliding photons of high energy, although, as mentioned earlier, such collisions naturally arise as a background 
in e + e~ collisions. It also reflects the difficulty in obtaining and presenting in a compact form complete analytical 
results, including bremsstrahlung, when massive fermions are present. The task would be even more formidable for 
reactions such as gg — ► QQ + X [[7|,|| . Our hope is that the approach of using high order series expansions to simplify, 
and clarify, such results will become more widespread. As well, our analytical integration of the single integral part 
is an essential part of a complete analytical integration, which is likely to be performed sometime in the near future. 
In the meantime, our parametrizations provide sufficient accuracy to be useful, as do our series expansions closer to 
threshold. 

The possibility of directly producing fermion pairs in 77 collisions has recently been realized at SLAC [^) , where a 
high energy (< 29.2 GeV) photon beam was collided on a low energy (2.35 eV) laser beam. Since the center of mass 
energy was insufficient to produce a pair (e + e~, in this case), multiple photon fusion was required; a different mech- 
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anism than that being considered here. The high en- 
ergy beam was produced via backscattering (of the same 
2.35 eV beam) off a 46.6 GeV electron beam and rep- 
resents a positive first step towards the construction of 
a higher energy 77 collider with both beams produced 
via backscattering. Of course, many technical difficulties 
arise in such a machine and these have been investigated 
(see M). 



II. GENERAL FORM AND DECOMPOSITION 
OF THE DIFFERENTIAL CROSS SECTION 

The process under consideration is 

7 (pi, Ax) + 7 (p 2l A 2 ) -> /(p 3 ) + /(p 4 ) + [V(k)], (1) 

where Ai, A 2 denote helicities and the pi, k denote mo- 
menta. / (= q, I) represents a fermion with mass m and 
V = g,j. The square brackets represent the fact that 
there may or may not be a gluon/photon in the final 
state. We have the following invariants, 

s=(pi+p 2 ) 2 , t = T -to 2 = (pi~p 3 ) 2 -m 2 , 
u = U — m 2 = (p2 — P3) 2 — m 2 (2) 

and 

s 2 = S2 — to 2 = (pi + P2 — P3) 2 — m=s + t + u. (3) 
Defining 



t 

s 



w ; 



s + t 



f3 = a/1 - 4to 2 /s, 
we may express 

t = -s(l - v), 
s 2 = sv(l - to), 

Now introduce 

a 2 e 4 f [7V c ] 
K ( a ) = 2tt ^— - 



l + (3' 



(4) 



-svw, 
s 



(i-n 



Z7T 



(5) 



(0) 



where the iV c , Cp factors are present only for / = quark 
and V — gluon, respectively and et is the fermion's frac- 
tional charge. Here 



ay 



a„ V 



Then 



a, V = 7 



da 

dvdw dvdw dvdw 
1 #(°> 



k(s) 



27r dudw 7r efociu> 



(7) 

(8) 
(9) 



The / functions are dimensionless functions of v and to, 
which allow us to parametrize our cross sections in an 
exact fashion, without dependence on ay. We use the 
normalization convention of jTl|| . Since, in that normal- 
ization, the contain an overall factor of n, we con- 
sistently present analytical results for f^/ir in order to 
cancel it. The unpolarized and polarized are given 
by 



/$p = ^[/ w (+.+) + / {0 (+,-)], 

£ ) i = J[/ (i) (+,+)-/ (i) (+,-)] ! 



(10) 



in the notation /W(Ai, A 2 ). Define 

j = l-<A 1 A 2 > 1 (11) 
where < Ai A 2 > is the average value of Ax A 2 . Then 

-f {l) (+-+) 

-r i> (+,-) (12) 



/ (i) (i) = 1± ^- A2 



1- < AiA 2 > (i) 



+ (+,+)] 



/«(+,+) 



(i) 



(13) 



so that 



.7 = 



- / w Cj) = / w (+,+) = / (i) (- -) 

- / (i) (i)=/ (i) (+:-) = / W (-,+) • (14) 

1 - / (i) (i) = /, w 



imp 



The LO term is given by 

1 



2-7T dvdw 



5(1 -w)< 



2m 2 /s 
v 2 (l — v) : 



■(l-2m 2 /s) (15) 



-J 



1 



2m 2 



v(l — v) 
S{1 - w) \J(l-v)z 

1 



-J 



v(l — v) 



1 



sv (1 — u) 



(16) 



2u(l - u) 



The last form shows explicitly the polynomial structure 
of the leading order differential cross section in terms of 
(3. This is somewhat misleading, however, as we shall 
see in the next section, since the phase space in v itself 
depends on (3. 

From |g] we see that df^ /dvdw has the form 

1 F fe („,«,)+ +Fs (v)5(l-w) (17) 



7r dudw 



= F h (v,w) + 



(l-w)+ 
F s (v,w) - F s (v,l) 



w 



2 



+F s (v,l)]n(l - wi)6(l - w) 
1 dfW 



7r dvdw 



1 d/W 1 d/W 
7r dvdw 7r dvdw ' 



where 



1U! 



l~ff 2 
4v(l - i>) 



and 



ld/W F s (v,w) 



7r dudiu (1 — w) + 



F s (v,w) - F s (v,l) 



1 — w 



(18) 



(19) 



(20) 



+F a (v,l)]n(l-w 1 )S(l-w) 



When integrating from w\ to 1, dp) /dvdw makes no 
contribution, otherwise it contributes. The function 
Fs(v)5(l — w) has the form 



F s (v)5(l -w)= F s (v, 1) In ( -j ) 5(1 - w) 



(21) 



2ir dvdw 
and F s (v,l) has the form 

1 



' ' ll<0) F 5 B ((3)+F s NB (v)S(l-w) 



F s (v,l) 



dp) 



2tt 3 v ~' dv 
Putting all these together yields 

1 dp) 



(22) 



7r dvdw 



F h {v,w) + 



F s (v,w)-±F S B({))^ 



1 — w 
ln(l — w\) + In 



4v 



1-/3 2 



2tt dvdw 

-&lfl + Fr{v)s{1 _ wh (23) 
27r dvdw 



Writing 



d-udu; \ dvdu; / 5 \ dudw ,/ ^ 



(24) 



where the subscript 5 denotes the part proportional to 
S(l — w), we have 



if HI) 

7r \ dvdw J s 



Fj((3) df^ 



2tt dvdw 



ln(l — wi) + In 



4v 



+ 



F s B (P)dp +F NB {v)S{1 _ w) (25) 



2tt dvdw 



and 

1 (dfW\ 
n \dvdwj Ng 



F h {v,w) + 



F s (v, W )-±F*(P)^ 



1 — w 



(26) 



The two simple F's are given by 
"1 + /3 2 




\nx + 1 



2mx - 2 



(27) 



Li 2 



-4/3 



(i - py 



In 2 a; 



-2 In a; In /3 + 2Li 2 (-a;) - 2Li 2 (x) - — 



(28) 



The other three F's are somewhat lengthy and will not be 
presented here as they can be directly inferred from the 
expressions given in O. F^ 8 is the contribution from 
the virtual diagrams which is not proportional to the 
Born term. Fh(v, w) is proportional to the first bracketed 
term in Eq. (30) of || and Fs(v, w) is proportional to the 
second bracketed term of that same equation. Both arise 
from gluon/photon brcmsstrahlung. 

It is standard || to divide the cross section (i.e. J*- 1 -*) 
into two parts. Firstly, there is the virtual plus soft part, 



v+s 



dp ] 



df. 



(i) 



dvdw dvdw dvdw ' 



(29) 



where dp) /dvdw denotes the virtual contribution and 
dp / dvdw is obtained by integrating the bremsstrahlung 
contribution to df^/dvdw over the region 

m 2 5 



1>W> Wi so ft = 1 



(30) 



then multiplying by <5(1 — w). We follow the defini- 
tion of the soft parameter, 5, given in [|| such that the 
gluon/photon radiated becomes arbitrarily soft by mak- 
ing 6 arbitrarily small. Since this takes into account all 
virtual corrections and soft radiation, the hard radiation 
may be taken into account by integrating df^ 1 ) /dvdw in 
the region w\ < w < wi. so [t- Since we never reach w = 1, 
F s (v, w)/(l — io) + = F s (v, w)/(l — w) in the hard radi- 
ation integration. We thus define df^ /dvdw as being 
df^ 1 ) /dvdw in the region w\ < w < u>i ;SO ft, where the 
S(l — w) terms do not contribute: 



dp 



dp 



(Wi < W < Wl, so ft). 



(31) 



It is not necessary to define dP /dvdw outside that re- 
gion since it is never evaluated there. 

The sum of the (integrated over some region) hard and 
soft contributions so defined is independent of S in the 
limit S — > and this method of separation is referred 
to as the phase space slicing method. As one might 

expect, there is a close relation between d/y+g/ 'dvdw 
and (df^ 1 ) / dvdw)$ as well as between dp) /dvdw and 
{dp) /dvdw) NS - We now give explicitly the necessary 
conversion terms. 
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It is straightforward to show that das /dvdw is ob- 
tained from the term proportional to doxo /dvdw in Eq. 
(30) of H by making the substitution 



In 6. 



(32) 



is 



From this, we infer that the conversion term necessary to 
transform our result into fyl s 

1 df(°) 



1 ^/sjonv 



7T dvdw 



2n dvdw 



In (5 — In 



ln(l — wi) 



4v 



1-ft 2 



The transformation is simply 

lffffs = 1 /#"»\ 
7r dvdw it \ dvdw J s 



1 df (1) 

+ "JS.conv 

7r dvdw 



(33) 



(34) 



To convert from (df^fdvdw)NS to df^ /dvdw we must 
add the following conversion term, 



1 dfi 11 1 /d/W 



7r dvdw 7r \ dwdw / NS 
1 



(1) 
conv 



7r V dvdw J 



N5 



7r dvdw 
1 Fj(ft)df^/dv 

2?7 1 — It) 



(35) 



Note that both soft and hard conversion terms are 
correctly determined by taking df^/dvdw into ac- 
count. Using (|32|), we reproduced the result of || for 
d/g] /dvdw and using (|34|), we reproduced the result of 

Q for dfy^'g np /dvdw. Further checks on df^ /dvdw will 
be discussed in the next section. 

The variables v,w are suitable for performing analyti- 
cal integration of the cross section (at least for the single 
integral part). They are not, on the other hand, suitable 
for performing series expansions of the integrated cross 
section about = 0. The reason is that, in these vari- 
ables, the integration limits depend on ft so that the series 
expansion of the integrated cross section does not follow 
straightforwardly from the series expansion of the differ- 
ential cross section, and we only have complete analytical 
results for the differential cross sections. Otherwise, we 
could just expand the final integrated result. The above 
will become clear in the following sections. This approach 
also allows for cross checking; when one first expands the 
differential cross section and then integrates, the result 
should coincide with that obtained by directly expanding 
the analytically integrated cross section. We will check 
this requirement for the single integral part, for which we 
do have analytical results. 

At this point, we introduce a new set of variables, r 
and w, suitable for performing series expansions in ft. 
They are defined through 

v=±(1+t0), w=l^c p (T)ft 2 (l-iu), (36) 



where 



We note 



1 - T 2 ft 2 



dfW _ /?3 C/3 ( T ) 



drdoj 



dvdw 



(37) 



(38) 



Because of the factor c^(r), df^/drdu) will never get a 
part proportional to 6(1 ± t). Defining 



Cu = Cf3(T)(l-Uj)(l + Tft), 

the invariants are given by 



(39) 



t 



\O--T0), 



U : 



'■(l + T/3-ft 2 C u ), 



S2 = ^C u ft 2 , 



T=--(l-2rft + ft 2 ), 



U 



^[l + 2r/3 + /? 2 (l-2c u )], 



(40) 



S 2 = -[l + ft 2 (2c u -l)}. 

In terms of r and u>, the LO term has the form 

1 rffW(j) ftS(l-w) 4 2 2 

2^ ~d^T = (1 - tW [2(1 -fi)-^ + r0) 

x(l + r 2 ft 2 -2ft 2 )] (41) 
= 5(1 - w)[(2 - j)ft + (2j + 4r 2 - 4jr 2 )ft 3 
+0(ft 5 )}. (42) 

We see explicitly that the j = 0, 1 differential cross sec- 
tions, in terms of these variables, vanish by order ft in 
the limit ft — > 0, while the j = 2 differential cross section 
is order ft 3 . 

The conversion term (B3I) becomes 



1 ^/s,conv 

7r drduj 



1 # (0) ^ 



2tt drdui 



Fsifi) 



In 5 — In 



2(l + r/3) 
1 -ft 2 



-\n[ft 2 cp(r)] 

and the conversion term ([55]) becomes 

1 Fi(ft)dfW/dr 



1 ^/h,coiiv 



7t drdco 



2tt 



1-LU 



(43) 



(44) 



III. ANALYTIC INTEGRATION OF THE DELTA 
FUNCTION PART 

The only complete analytical results for the differential 
cross sections were presented in Analytical results 
for the virtual+soft part were presented in || for the 
unpolarized case and in [|| for the polarized case (where 
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the virtual and soft parts are given separately, in terms 
of various functions). Still, not even the virtual+soft 
part has previously been integrated (over fermion angle) 
analytically. In this section, we present such an analytical 
integration. We were not able to integrate the non delta 
function (or hard) part analytically, in a straightforward 
fashion, and reserve that for future work. 

The integrated cross section (or f( l >) is obtained via 



(0 



/ 



where 



dv / dw 



df 



(0 



dvdw 



1 f 1 d/W 
dr / du}—— , 
i J drdu 



Vi = -(l-(3), 



v 2 = -(! + /?). 



(45) 



(46) 



Let 83 be the angle between p 3 and p\ in the 77 cm. 
Then #3 is given by 



1 — v — vw 



cos O3 



-V + Vw) 2 + (3 2 - 1 

2r - /3c„ 



^A-c u (4-(3 2 c u ) 



Thus, 



COS 6*3 



1 - 2v 



for w = u> = 1. 



(47) 
(48) 

(49) 



We see that for (3 — > 0, cos #3 varies rapidly with while 
it is simply equal to r. This is why the phase space in 
v becomes vanishingly small by order [3. Similarly, from 
([|5|) and ([L9|), or (36), we see that the w phase space is 
order /3 2 . Thus, the double integration over v and w is 
order /3 3 , in accord with (p8[). 

The integration of (|l|) or @ is rather straightfor- 
ward, yielding the LO term, 



2tt 



2/3(1 + /3 2 ) - 6/3 j - (1 - /3 4 + 2j) lax. (50) 



Since 



00 a2k+l o3 

**=-*Y,wri = - 2p - 2 T-'"> (51) 



we have 

/ (0) 0') 

2tt 



2(2-^+^/33 



(52) 



fc=2 



2fc-3 2fc+l 



so that / (0) (0,1) are order (3 and / (0) (2) is order /3 3 . 
Also, we see that only f^°'(0) is finite in the limit j3 — > 1 
and it approaches 



/ ( °)(0)^8tt, for/3^1 



(53) 



This is because the 1 — /3 4 term in (^l]) keeps the j = 
channel finite. For j — 2, the cross section vanishes for 
exactly t = ±1, as required by angular momentum con- 
servation along the 77 axis, but for /3 — ► 1 the part pro- 
portional to j goes like (1 + t 2 )/(1 — r 2 ) as soon as we 
move away from exactly r = ±1 and is hence not inte- 
grably finite at j3 = 1. In order that the j = cross sec- 
tion be nonvanishing for r = ±1, where its maximum lies, 
the / and / must have opposite spins by angular momen- 
tum conservation, leading to m J s (~ 1 — /3 2 ) suppression 
in the numerator. The fact that the LO j = cross sec- 
tion continues to be 1 — (3 2 suppressed for r ^ ±1 follows 
from symmetry arguments [Q . This exactly compensates 
the t-channel singularity in the propagator, leading to a 
finite / (0) (0) for /3 -> 1. Of course, for t ^ ±1, the j = 
differential cross section will vanish like (1 — /3 2 )/(l — r 2 ) 2 , 
making it unobservable in LO, for /3 — > 1. So, had we 
taken the limit /3 — > 1 from the beginning, the j = 
cross section would have vanished identically. Hence the 
nonzero /^(O) in the /3 — ► 1 limit is a remnant of using 
the fermion mass as a "regulator" . 

Near threshold, the 1 — /3 2 suppression of the j = 
channel will not be significant, hence the major con- 
straint will come from angular momentum conservation 
in the forward and backward directions which will lead 
to suppression of the j — 2 cross section there. The 
j = cross section reaches its maximum in those config- 
urations, however. Thus, we can clearly understand the 
feature of the numerical results for top quark production 
in Q which show that imposing angular cuts in the di- 
rection of the beam pipe has a greater effect on the j = 
channel than on the j — 2 channel. 

We denote the single and double integral contributions 
to /« by 



(i) 

si/ di 



fir I dw ( ^—r-\ 

\dvdwJ s/NS 



Vl 

dT [ 

-1 Jo 



duJ \dTdu) s/NS ' 



(54) 
(55) 



We performed the single integration using (|34]), as op- 
posed to (|5^). It turned out to be quite lengthy and in- 
volved. We did not check to see whether using (^5|) sim- 
plifies the calculation. This question is probably more 
relevant to the double integration, however. Our final 
(simplified) result is 



f 1 \n(x) 
— = ai n 2 /6 + a 2 U 2 (x) + a 3 Li 2 (-a:) + — ^ 

7T (j 



a4Li 2 (x) 



+a 5 Li 2 (-x)} + a 6 ln[(l + /3)/2][tt 2 /6 + 2Li 2 (-a?)] 
+a 7 ln(x)7r 2 /6 + a 8 {-31n 2 [(3 + /3 2 )/4] ln(x) 
+51n[(l -/3 2 )/4] ln[(3 + /3 2 )/4] ln(a;) 
+2Li 3 [-2z/(l + /3)] - 2Li 3 [-2/( a; (l - /3))] 
+2Li 2 [(l - /3) 2 /(3 + I3 2 )} ln[2/(.x(l - /3))] 
-2Li 2 [(l + /3) 2 /(3 + /3 2 )] m[2x/(l + 0)]} 
+<z 9 {Li 3 [(l + /3)/2] - Li 3 [(l - /3)/2]} (56) 
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+a w ln(ir) + an ln 3 (x)//3 + a%2 ln 2 (x)//3 

+ai 3 ln 2 [(l + /3)/2] ln(x) + a u /3 

+ [a 15 + a 16 ln 2 (x)//? 2 ]/?ln(/?) 

+ [ai7/3 + a ls ln(x) + aig ln 2 (cc)//3] ln[(l + 



/3)/2] 



The coefficients a,i(j) are given in Appendix A and are 
rational polynomials in /3, finite as (3 — > 0. Terms of the 
form F(/3) — F{—(3) will vanish for /? = 0. Noting that 



Li 2 (l) = T , 



Li 2 (-1) = -^, 



(57) 



we see that only the terms proportional to a\, . . . , 05 
may contribute at threshold, which is the case for j = 0. 
Indeed, one finds the correct threshold correction from 
those terms alone. The relevant series expansions will be 
given in the next section. As we shall see in Section V, 
the double integral series expansion starts at order /3 3 . 

Two independent determinations of (pq ) were made us- 
ing Mathematica fl§ and REDUCE |sf That software 
could not evaluate certain integrals which can be found 
in |u| . It was verified that the analytically integrated re- 
sult agreed numerically with the numerically integrated 
result. In the next section, we will show how one can use 
the series expansion as a very solid check as well. 

Perhaps the most convincing check of (pq) and the 
analytical result for df^/dvdw (or da^ho! dvdw) ob- 
tained in p| is the excellent numerical agreement with 
tabulated results for /W existing in the literature. The 
only existing analytical results, aside from those in 

are the expressions for df^l 3 /dvdw (i.e. /dvdw)g), 

dfg/dvdw given in Q for the unpolarized case (us- 
ing dimensional regularization), with which we agree 
exactly, and similar expressions for the polarized case 
in H (obtained using a gluon energy cut and a small 
gluon mass as infrared regulator). The latter are not 
quite in a form suitable for direct analytical comparison. 
There have been no other analytical results presented 
for (d/W / 'dvdw)N8 in the polarized or unpolarized cases. 
Hence we must perform the above mentioned numerical 
checks. 
Define 



2m 



1 



P = y/1- 1/z 2 



(58) 



In Table | we give numerically computed values for / U n P) 
/poi' /^(+> +)> ~) as weu as the specific contri- 

butions from all the fty and fff to the corresponding 
/W, for various values of 1.2 < z < 20. The result at 
z = 1 is given exactly by the series expansions presented 
in the next section. We also indicate the number of sig- 
nificant figures, n.s., following the decimal point, in 
(and /<')). 

We find it useful to describe how the values in Ta- 
ble I were obtained in order that one may see clearly 



The val- 
9) which 



which numbers have been rounded and how. 
ues for / s ( ^ unp , /f* pol were obtained using 
permits arbitrary precision, using a package like Math 
The values of f£\+,+) and f£\+,- 



ematica. 



were obtained adding/subtracting the values of f siu 



/iiiol so obtaincd - Thc /il.lnp' 



f (i) 

Id 



(1) 

52,u: 

idi,po\ were obtained 
by numerical integration using (Q) for thc NS part. 
From these, f^(+,+) and fjil\+,—) were obtained by 
adding/subtracting. Finally, / (1) (+,+), / (1) (+>~) 

were 

obtained by adding the corresponding f£' 



than adding/subtracting /„n P , /po|); similarly for /, 



fdi (rather 



c(l) 

/ unp , 



/poi- ^ e not cnec k to see if (^5|) leads to any reduc- 
tion in computational time, for the precision obtained. 
The general trend is that one needs more integration 
points as one goes to higher z. 

The next issue is, of course, how well these values com- 
pare with other tabulated values for /W . Two other such 
tables exist at present. The original one of [jllj gave /unp 
for z = 2, 3, 4, 5, 10; the value at z = 1 being numerically 
equal to the known threshold result, as given in the next 
section. Their numerical values were obtained using the 
/v+s nP gi yen in H), added numerically to unp , deter- 
mined there using the same methodology as |] , which is 
equivalent to our method. We find numerica agreement 
with pl[ | to within the precision of those values, which is 
roughly at the order of one part in 10, 000 or better. This 
can only be achieved with correct analytical results. Our 
calculation of is identical in method (same integrals 

and structure) to that of / un p (at the differential and in- 
tegrated level), the only difference arising from different 
traces due to the contraction with a polarized photonic 
tensor rather than an unpolarized one. As two indepen- 
dent determinations of these traces were performed, there 
is little room for any error in Z^,]. Fortunately, we may 

directly check this assertion since the values of/W(+,+) 
and /W (+, -) for z = 2, 3, 4, 5, 10, 20, 50 were tabulated 
in ||. There, Monte Carlo methods were used, leading to 
accuracy at the level of better than 1% in regions where 
the are sizable, but apparently not better than ±0.2 
or so in absolute error. This absolute error is noticeable 
only for / W (+, +) and only for z — 2, 3, where f^'(+,+) 
is small. To within the above accuracy, we are in good 



agreement with ||. Since / (1) (+,-) = Amp- / po } and 
since we have precision agreement with fll| for /u np and 



with 



for /«(+,- 

' po 



-), we conclude that our analytical 

results for df^/dvdw of Q have been verified. In light 
of the above, Table | is seen to be the most complete and 
precise such table at present . 

We may convert from ( fitf ) to /y+g by adding the fol- 
lowing conversion term, 



f (i) 
Jv+s 



f(l) 



J S.conv 



(59) 



G 



(1) 



In 



In S 



2tt 



+(3(1 + 1 - 4j) + 2/3(1 + 1 - 3j) In 



1 + /3 



2/3 2 

~ + (3 - 2(3 2 + (3 3 - tL. + j(4 - 3/3 + /3 2 ) 



•In a; + 2 In x In 



/? 4 
2 

1 + /3 



In a; 

2Li 2 (x) - Li 2 (-a;) 



+ - 



(l-/3 4 + 2j) , 



where /s C o nv follows from integrating df^ onv /dv, given 



(i) 



in (33), or from integrating df^l onv /d,T, given in (|43| 



IV. SERIES EXPANSION OF THE DELTA 
FUNCTION PART 

Besides providing a useful check of the analytical inte- 
gration of the previous section, there are many reasons 
why it is useful and instructive to series expand the dif- 
ferential and integrated cross sections about (3 — 0. In 
the absence of complete analytically integrated results, 
only a series expansion about (3 = can be used to make 
(very) high precision predictions in the (3 ~ region. 
One also sees the structure of the cross section in a way 
that cannot be inferred from the non-expanded analytical 
results, which are somewhat complicated. From a practi- 
cal viewpoint, having "simple" series expansions for the 
differential cross sections allows one to do complete nu- 
merical studies in the region not too far above threshold 
rather easily. This is because the resulting expansions 
only involve simple polynomials and simple logarithms. 
We will address the issue of the region of validity of the 
expansions as well. 

The other issue is that of resummation. There are large 
correction terms at threshold which can be resummed. 
Having a series expansion of high enough order to be of 
practical use allows one to explicitly perform resumma- 
tions up to some order in [3 while leaving the higher order 
terms the same. The net result would be an equally sim- 
ple series, improved via resummation so as to allow one 
to go closer to threshold. This is beyond the scope of 



this paper as are other very near threshold effects. Suffice 
it to say that having the threshold series expansion will 
facilitate these studies for those interested. 

Throughout, we will expand up to order (3 10 (includ- 
ing (3 11 h\(3 terms). The expansion which exists in the 

literature (see [|ll|) is only for /unp and only goes to or- 
der (3. Going to order (3 10 may seem excessive at first, 
but we found it to be a good stopping point for several 
reasons. Considerable structure arises beyond order (3 
which allows us to see the general, all-orders in (3, struc- 
ture of the various series. Also, one gains little in terms 
of precision by going to even higher orders in (3, with- 
out including several more terms. Then, the series would 
start to become lengthy and cumbersome, reducing the 
advantage over the analytical result in terms of ease of 
use. For certain series, going much beyond (3 W would 
take a very large amount of computer memory and run- 
time, not justifying the extra effort, as going to order (3 
was a considerable task in itself. Finally, by going to such 
a high order, we may stringently check the analytically 
integrated single integral result of the previous section as 
will be described below. 

We find that df^ /drdoj may be expanded in the gen- 
eral form 



df~W 
drduj 



(60) 



i=0 j=0 

Therefore J*- 1 ) may be expanded as 

oc 1 

/ (1) = / (1) = EE^ lnJ ^ 

i=0 j=0 

where the rfy are given by 



dr / dui Cij(r,uj). 



(61) 



t-i Jo 



(62) 



With the variables v and vo, the integration limits depend 
on /3, hence the above arguments do not hold. So, one 
sees clearly the necessity of the change of variables. 

We convert from (dfW/dvdw) s to (dfW/dTdui) s using 
(pq), which modifies the overall factor via 5(1 — w) — > 
/3 3 c^(r)<5(l -w)/2 = (36(1 - lu)/2. Then, the results for 
the series expansions of (df^ /drdu>)s are, for j = 0, 
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1 ( df {1 \+,+) 
7r V drdui 



6(1 - w)|2tt 2 + /?(-20 + 7T 2 ) + 2/?V(l + 2r 2 ) + /3 3 (-9tt 2 (2 - 7r 2 ) - 4(22 + 45r 2 ) 

+ 96{(1 - 9r 2 ) ln(2) + ln[4/3 2 (l - r 2 )]}) /9 + 2/3 4 [-3tt 2 + 16r + 3tt 2 t 2 (2 + 3r 2 )]/3 

+ /3 5 (225tt 2 (1 - 14t 2 + 17t 4 ) + 4[1394 - 25(464r 2 - 381r 4 )] + 480{(-34 + 205r 2 - 245r 4 ) ln(2) 

+ 2(1 + 5t 2 ) ln[4/3 2 (l - r 2 )]})/225 + 2/3 6 [16r(6 + 35r 2 )/45 - tt 2 (1 + r 2 )(l + r 2 - 4r 4 )] 



/3 7 
32 



7T 2 r 2 (7 - 34r 2 + 31t 4 ) + 4{-72244 + 49r 2 [35068 - 5r 2 (18680 - 12373r 2 )]}/11025 



>({227 + 7r 2 [-473 + r 2 (1315 - 973r 2 )]} ln(2) + [-26 + 7r 2 (4 + 15t 2 )] 
x ln[4/3 2 (l - r 2 )]) /105] + 2/3 8 {7r 2 r 2 (-l - r 2 )(2 + r 2 - 5r 4 ) 
+ 32r[-39 + 7t 2 (7 + 29r 2 )]/315} + /3 9 3157r 2 r 4 [17 + r 2 (-62 + 49t 2 )] 
+ 4(830486 + 9r 2 {-2777328 + 7r 2 [2146844 + 65r 2 (-55762 + 29125r 2 )]}) /315 



+ 32^-376 + 3r 2 {4899 + r 2 [-26579 + r 2 (46207 - 25243r 2 )]}j ln(2) + 64{-ll + 3r 2 [-26 
+ 7r 2 (3 + 10r 2 )]} ln[4/3 2 (l - t 2 )]] /315 + 2/3 10 (^ 2 r 4 (-l - r 2 )(3 + r 2 - 6r 4 ) 
+ 32r{-55 + r 2 [-455 + r 2 (406 + 1455r 2 )]}/1575) 

+ 64/3 11 ln(2/3) (-122 + llr 2 {-44 + 3t 2 [-78 + 7r 2 (8 + 25t 2 )]}) /3465 j 
and, for j = 2, 

ir \ drdu) 



2(1 - t 2 )S(1 - w)|2/3 2 ^ 2 - 16/3 6 + /3 4 vr 2 (l + 5r 2 ) 

+ 4/3 5 {14 - 159r 2 - 6(7 - 19r 2 ) ln(2) 

+ 241n[4/3 2 (l - t 2 )]}/9 - /3 6 {-32r/3 + tt 2 [1 

- t 2 (3 + 8t 2 )]} - 2/3 7 (ll02 - 21355t 2 + 37995t 4 
+ 120{(-19 + 273r 2 - 382t 4 ) ln(2) 

+ (1 - 25r 2 ) ln[4/3 2 (l - t 2 )]}) /225 

- /? 8 r[16(3 - 85t 2 )/45 + tt 2 t(2 - 5t 2 - llr 4 )] 

+ 2/3 9 {129611 - 805r 2 [3721 - 2t 2 (6781 - 6092t 2 )] 
+ 840[-219 + t 2 (5172 - 18803t 2 + 16134t 4 )] ln(2) 
+ 3360(1 + 70t 4 ) ln[4/3 2 (l - t 2 )]}/11025 
+ /3 10 r[7r 2 r 3 (-3 + 7T 2 + 14r 4 ) 
+ 16(12 - 7t 2 + 1057r 4 )/315] 

+ 32/3 11 ln(2/3)[13 + 51r 2 + 21r 4 (l + 55r 2 )]/315 



(63) 



(64) 



The remaining terms are of order [3 11 . Here, and throughout, we group terms proportional to ln(2/3) rather than ln/3 
as in (|60|), (61) for purposes of compactness. Of course, the choice is quite arbitrary. 

We notice that the cross section is isotropic up to order /3; the angular (t) dependence enters only at order (3 2 . The 
LO term, on the other hand, was isotropic up to order 1 . We also see that the step function threshold behaviour 
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arises entirely from the j : = channel, at the level of 
the differential cross section, since the j = 2 channel 
starts at order 1 . From the 1 — r 2 overall factor, and 



using (^9() , we see that the delta function contribution to 
the j = 2 cross section vanishes at cos #3 = ±1 as did 
the LO cross section (^). This vanishing is not obvious 
from the exact analytical expressions, but simply reflects 
angular momentum conservation along the 77 axis when 
u> = 1 (2 — > 2 kinematics). The j = channel, on the 
other hand, becomes infinite (but integrably finite) for 
cos #3 = ±f due to the l n(l — t 2 ) terms. 

The expansions (|6^), ( |64| ) have rather simple structure 
in that, aside from the ln(l — r 2 ) terms, the Cy (r) are 
simply polynomial in r. This amounts to considerable 
simplification and reduction in computational time rel- 
ative to the exact expressions, especially after the non 
delta function part is added, where the simplification is 
even greater as we shall see in the next section. Adding 
the conversion term (^|) gives ^dfyK/drduj. 

Two independent calculations of (|65|), ( |64| ) were per- 
formed using Mathematica and REDUCE. The expan- 
sions were also checked numerically by subtracting them 
from the exact expressions. The difference was checked 
to be of order This is most straightforwardly done 
by taking rather small (3. 

Assuming we are working at (3 where the series are suf- 
ficiently accurate, one could easily analytically integrate 
(|63|), ( p4| ) over a region of r {cos 6 3) relevant to some 
experiment, if desired, and implement angular cuts ana- 
lytically. After a suitable change of variables, the same 
could be done for the hard radiation part, either ana- 
lytically or numerically. Cuts on additional observables 
may be made by subtracting off the unwanted configura- 
tions using the squared amplitudes given in §| and Monte 
Carlo integration, for instance. Here, we simply present 
the total integrated results. 

For the j = channel, we find 



f 



-/£>(+,+) = 

TT 

2I 2tt 2 - (20 - tt 2 )/5 + 10/3tt 2 /3 2 + /3 3 /3[- 340/3 + tt 2 



+641n(2/3)] + 8/1 5tt 2 /3 4 + 4/ 15/? 5 [-6343/45 - tt 2 
-32 ln(2) + 256/3 ln(2/3)] - 104/105tt 2 /3 6 
+4/105/3 7 [-39163/315 - tt 2 + 208/3 ln(2/?)] 
-88/3157T 2 /? 8 + 4/9/? 9 {1/7[-tt 2 /5 - 64/3 ln(2)] 
+l/25[1281n(2/3) - 43903/315]} - 488/3465tt 2 /3 10 

+103232/51975/3 n ln(2/3) 



(65) 



and for j = 2, 



1 



(66) 



+ ln(2/3)/3] + tt 2 /7/3 6 + 6/5/3 7 [-3947/945 + 16/7 ln(2) 
+32/9 ln(2/3)] + 29/105tt 2 /3 8 + 4/15/3 9 [-823/45 + 8 ln(2) 

+16 m(2/3)] + 289/11557r 2 /3 10 + 256/63/3 11 ln(2/3) 

The results are indeed quite simple. We may obtain 
~/v+s ^ adding the conversion term (p39|). 

The strongest check comes from the fact that the ex- 
pansions (|65|), ( |66| ) which come from integrating (j6^), 
(|64| ) agree exactly with the expression obtained by ex- 
panding the analytically integrated result d56| ) directly. 
In this way, we simultaneously check all the above men- 
tioned expressions, including our analytical integration 
(f56|). The expansions (|65|), ( |66| ) were also checked nu- 
merically by subtracting them from the exact expression 
d56| ) and verifying that the difference was order /3 11 , as 
was done for (|63|), (B4l). 



V. SERIES EXPANSION OF THE NON DELTA 
FUNCTION PART 

Perhaps the most remarkable result of the series expan- 
sion is the simplification of the non delta function part, 
whose original form is the most lengthy part of the exact 
result, involving complicated logarithms, etc. . . Although 
the intermediate expressions were very lengthy and con- 
siderable computational time was required, a large degree 
of cancellation resulted in the following simple series. We 
convert from {df^'/dvdw)NS to {df^ 1 '/drdu))Ns by mul- 
tiplying by P 3 cp{r)/2. Then, for the j = channel we 
find 



1 / rf/W(+,+) 
tt \ drdco 



8/3(1 - r 2 )i 4/3 3 + 4r/3 4 - l/5/3 5 [27 - 43w 



16/3<^ 7r 2 /r - 8/3 3 + tt p + 32/3° [-289/720 + ln(2)/5 



-t 2 (63 - 43w)] - 2/5t/3 6 [15 - 43w - r 2 (53 - 43w)] 
+ l/35/3 7 [2(-87 - 108w + 179cj 2 ) - r 2 (1133 - 2353w 
+716u 2 ) + t 4 (1499 - 2137w + 358w 2 )] 
+8/35r/3 8 [70 - 259w + 199w 2 - 2r 2 (133 - 416w 
+ 199w 2 )+r 4 (339- 573w + 199w 2 )] (67) 
-2/105/3 9 [338 - 1363w + 2222w 2 - 1189w 3 
+r 2 (-1639 + 11135w - 13399w 2 + 3567w 3 ) 
-r 4 (-8663 + 27244w - 20132w 2 + 3567w 3 ) 
+r 6 (-7606 + 17472cj - 8955w 2 + 1189w 3 )] 
+2/105r/3 10 [-462 + 4317w - 9396w 2 + 5221w 3 
+3r 2 (2529 - 12399w + 15371w 2 - 522kj 3 ) 
+3t 4 (-6022 + 23659cj - 21346w 2 + 5221w 3 ) 

+r 6 (13897 - 38097^ + 27321w 2 - 522W 3 )] j 
and, for j = 2, 
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1 ( dfM(+, 

7r V drdu) 



N5 

-8/3(1 - r 2 )|4/3 5 [l + uj + t 2 (1 - u)] - 8t/? b (1 - w) 

x (1 - r 2 ) + 2/5/3 7 [3(l - 6w + 7w 2 ) - 2t 2 (14 - 31w 
+21w 2 ) + r 4 (73 - 44w + 21w 2 )] - 2/5r(l - r 2 )/? 8 
x [-53 + 30w - 35w 2 + r 2 (169 - 154w + 35w 2 )] 
-l/35/? 9 [152 - 226w + 385w 2 - 415w 3 
+r 2 (-2688 + 3890w - 3451w 2 + 1245w 3 ) 
-r 4 (-7192 + 8782w - 5747w 2 + 1245w 3 ) 
+r 6 (-5792 + 5118cj - 2681w 2 + 415w 3 )] 
+2/35(1 - t 2 )t/3 10 [-420 + 1035w - 1475w 2 + 904w 3 
+2r 2 (2091 - 2569w + 2215w 2 - 904w 3 ) 

+r 4 (-6390 + 6971w - 2955w 2 + 904w 3 )] 1 . (68) 



We notice the absence of any logarithms, including 
powers of In [3. The structure is fairly predictable as 
well. We see that the series begin at order /3 3 and /3 5 
respectively, so that their effect will be negligible very 
near to threshold. On the other hand, the large coeffi- 
cients imply that they soon become noticeable for small 
P. We may obtain Id^/drdu by adding to @, 
the conversion term (|44|) . 

Two independent determinations of (|67|), ( |6"8| ) were 
performed using Mathematica and REDUCE. These ex- 
pressions were also checked numerically analogously to 
the delta function part of the differential cross section. 

The integration of (|67]), ( |68| ) over r, ui is straightfor- 
ward and we obtain 



7T 



-128 „ s 448 n * 34624 

B 6 B 5 H 

9 225 2205 



42368 



3675 



/3 9 , 



(69) 



I . Wfj _ _^ _ ^1024^5 _ 2816^7 _ 134656^ 



(+' ) - 45 



525 



19845 



B J . (70) 



These are remarkably simple results, which suggest that 
the exact integrated result for ffj^ is not too complicated. 
We notice the vanishing of the coefficients of the even 
powers of f3. This follows from the antisymmetry in r of 
the corresponding terms in the differential cross section. 
To convert to fjp , we must integrate the conversion term 
( |35| ) (or (pi)), over v and between w\ < w < iUi )S oft- This 
yields 



/h 



f 0) 

J di , ■'H,conv 



7T 

t 0) 



/•(I) 

J S,conv 



where /i 1 ' is given in (59). This verifies the cancella- 



tion of the S dependence of + in the limit 5 — ► 0. 



Implicitly we were working in this limit since we inte- 
grated the non delta function part over all r, u>. 

We checked ([39]), ([70]) against the numerically inte- 
grated result and again found the difference was order 
(3 , In the next section we will tabulate the numerical 
errors on the series expansions for /W (total) for vari- 
ous values of B, relative to the numerical result obtained 
from the exact expressions. 



VI. TOTAL SERIES RESULTS AND 
NUMERICAL PARAMETRIZATIONS 

We are now in a position to study the total cross sec- 
tion, by combining the results of the previous sections. 
Adding (|65| ) and ( p9| ) gives the series for the j = total 
cross section 



1 



/ (1) (+,+) 



2|2tt 2 - (20 - ir 2 )8 + 10/37T 2 /? 2 + /3 3 /3[-404/3 + tt 2 

+641n(2/3)] + 8/15tt 2 /3 4 + 4/15/3 5 [-6511/45 - tt 2 
-321n(2) + 256/3 ln(2/3)] - 104/105tt 2 /? 6 
+4/105/3 7 [25757/315 - tt 2 + 208/3 ln(2/3)] - 88/315tt 2 /3 8 
+4/9/3 9 {l/7[-7r 2 /5- 64/3 ln(2)] + 1/25[128 ln(2/3) 
+407639/2205]} - 488/3465tt 2 /3 10 



-103232/51975/3 n ln(2/3) 



(72) 



and adding (p6|), ( |70| ) gives the series for the j = 2 total 
cross section 

1 / (1) (+,-) = 

7T 

16/3|tt 2 /3 2 - 8/3 3 + 7T 2 /? 4 + 32/3 5 [-77/144 + ln(2)/5 

+ ln(2/3)/3] +tt 2 /7/3 6 + 6/5/3 7 [-677/135+ 16/71n(2) 
+32/9 ln(2/3)] + 29/105tt 2 /3 8 + 4/15/? 9 [-16949/735 
+8 ln(2) + 16 ln(2/3)] + 289/11557T 2 /? 10 



+256/63/? 11 ln(2/3) 



(73) 



Such simple expressions indeed make numerical studies 
not too far above threshold rather straightforward. We 
can get an idea of how well these series work for typical B 
by comparing with numerically calculated values of f^. 
In Table Q we present the fractional error on the se- 

(71) ries for / (1) (+, +), / (1) (+, -), fun P relative to the result 
obtained using numerical integration, for various values 
of z in the region 1.05 < z < 1.4. For z < 1.05, the 
series expansions are more accurate than the numerical 
results. At z = 1.05, the errors are at the 10~ 7 — 10~ 6 
level. For z = 1.2 they are at the 10~ 4 — 10 -3 level and 
for z = 1.4 they are at the 10~ 3 — 10 -2 level. The errors 
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on f^(+,+) are at the lower end, while the errors on 
/^(+, — ) are at the higher end and those for /iip lie 
in between. This is good because, as we shall see in the 
next section, in determining ot s via top quark production 
at a 77 collider, it is the j = and unpolarized channels 
which are of interest, the j = channel being the most 
interesting one. With precision of better than one per- 
cent for z < 1.4, we have sufficient accuracy to use the 
series expansions (differential in particular) to perform 
easy numerical studies relevant to top quark production 
at a 77 collider of ^/s < 500 GeV. As we shall see, for the 
a s determination, going to much higher energies is not 
useful since the determination is best done near z = 1.2 
{y/s ~ 420 GeV). 

It is also useful to be able to parametrize /W to good 
accuracy for larger (3, relevant for bottom and charm 
quark production at intermediate energies or top quark 
production at very high energies. This was done by fit- 
ting numerically computed values of f^. We divide 
the parametrizations into 3 regions: a low energy re- 
gion (1 < z < 1.5 or < (3 < 0.7454), an intermediate 
energy region (1.5 < z < 5) and a high energy region 
(5 < z < 20). We will denote the corresponding /W as 
fW' le , f^' le and f^' he , respectively. 

The various forms for the parametrizations are 



27r 2 -(20-7r 2 )/3+^/3 2 



64 o 
+ y/3 3 In/3 



i=3 



/(!).-(+,. 
/( !Me (+) . 



= ^ Cl (z-i.sr 

4 

= 5>(*-5) j 



(74) 



parametrizations to include resummation effects without 
changing the higher order coefficients. Here, we simply 
present the one-loop corrections. 

In the intermediate energy region, f^' le (+, — ) is accu- 
rate to < 0.1%. / (1) ' ie (+, +) is accurate to < 1%, except 
very near /W' ie (+, +) = 0, which occurs for z ~ 2.15, 
3.15. There, the absolute errors remain small, but of 
course the fractional error is larger. In the high en- 
ergy region, 

/(i),fce( +) _) i s accurate to < 0.05%, while 
f (1) ' he (+,+) is accurate to < 0.5%. The above errors 
are rather conservative and one can not distinguish the 
parametrizations from the exact results for practical pur- 
poses. 

The (exact) plots of f' l \ in the three energy 

ranges are given in Figures [l] [|. In Fig. |l] we plot 
f( ' in the low energy region versus (3. [3 is more su- 
tiable than z in this region since the threshold region 
becomes compressed and /W varies quite rapidly with 
z right at threshold. Fig. [I] (a) highlights the fact that 
/W is most naturally decomposed into f^(+,+) and 
/^(+, — ) since f^{+,+) is monotonically decreasing 
in the threshold region while — ) is monotonically 

increasing, /unp on the other hand exhibits a rather sud- 
den dip and peak which seems unnatural, until broken 
down into (+, +) and (+,— ) components. From Fig. [j] 
(b) we see that is monotonically increasing in this 
region for both helicity states. 

In Fig. |2| we plot f' l \ versus z in the interme- 
diate energy region. We notice that the /W cross just 
before z = 1.5 and the cross just before z = 2. 
f^(+,+) becomes negative and reaches a minimum 
near z = 2.6 and then monotonically increases, as does 
throughout. / (0) (+,-) 

continues to grow 
while / (0) (+,+) levels off in accord with (|53|). In Fig. 
| we plot /W, in the high energy region. The be- 
haviour remains unchanged. 



and 



/(!).'«(+,-) 



16tt 



f^( +) -) = J2ci(z -1.5)' 

i=0 

/«' fte (+,-) = 5>0*-5)\ 
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[^/3 2 -8/3 3 + 7r 2 /3 4 ] 



i=5 



(75) 



The Ci are given in Appendix B. In the low energy re- 
gion, where high accuracy is required, the parametriza- 
tions are accurate to < 0.01%, with the errors being the 
largest near the higher end of the region. The leading 
terms, given analytically, guarantee the correct thresh- 
old behaviour as they are just those in the exact series 
expansion. As mentioned earlier in connection with the 
series expansions, one can explicitly perform resumma- 
tions on those terms. Thus one could modify the above 



VII. PRECISION a s DETERMINATION FROM 
TOP-QUARK PRODUCTION 

A high energy 77 collider can be used as a "factory" 
for many interesting particles: Higgs bosons, bosons, 
top quarks etc... The beam polarization we be useful in 
producing Higgs bosons and reducing QQ backgrounds. 
More specifically, the j = channel will be of interest. 
This channel also turns out to be the channel of interest 
when trying to determine a s via top quark production, 
making it complementary to the Higgs studies. The rea- 
son is that the cross section, and QCD corrections, are 
enhanced in this channel, thereby improving the statis- 
tics and the determination of a s , to which the cross sec- 
tion will be quite sensitive. The process 77 — > tt + X 
is more powerful than e + e~ — > QQ + X in determining 
a 8 because the QCD corrections are quite small in the 
latter, thus requiring an unreasonably large number of 
events for high precision; the corrections are suppressed 
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by a s /ir ~ 4%, relative to the Born term. In 77 — > tf+Jf, 
we can "pick" our QCD correction by choosing the ap- 
propriate beam energy. Of course, as one gets too close 
to threshold, the perturbation series cannot be trusted, 
for reasons we will discuss below. Hence there are limi- 
tations. 

To best illustrate the above idea, in Fig. || wc have 
plotted the 77 — > it + X cross section at LO and NLO, 
in the region 1 < z < 1.4, for the various helicity states. 
We took N f = 5 , m t = 174 GeV and used A = 230 MeV 
in the two- loop expression for a s , evaluated at [i 2 = s. 
One could also use Nf — 6, but since we are not far above 
threshold it is simpler to use Nf = 5 for evolution from 



/i 2 = M| to /1 2 = s. We notice that the j = cross 



section is the largest, as are its QCD corrections, in this 
region. The region z ~ 1.2 is nice in that the j = cross 
section is near its maximum and the QCD corrections 
are sizable (~ 20% of the total cross section), yet not so 
large that the perturbative expansion is unreliable. As 
one gets closer to threshold, other higher order effects, 
nonperturbative effects and top width effects may also 
become important. For these, and other reasons to be 
considered below, we will suggest z = 1.2 as being the 
optimal region for extracting a s and we will give a rough 
estimate of how precisely a s may be determined there. 
As well, we suggest the j = channel as being the most 
powerful. 

Firstly, we note that z = 1.2 corresponds to ^/s 77 ~ 
420 GeV, for top quark production. This energy should 
be accessible at a y/s e + e - > 500 GeV NLC. A typical 77 
luminosity assumed is 20 fb _1 . Since er ~ 1.4 pb, this 
corresponds to roughly 28,000 it events. Since the QCD 
correction is ~ 20% of the total cross section, this trans- 
lates to Aa s /a s ~ 3%, statistically. With a luminosity 
increase and, possibly, extended running, one could en- 
vision going to the percent level or better. 

The above analysis was purely based on statistics and 
one-loop QCD corrections. Therefore, we will briefly dis- 
cuss various theoretical systematic uncertainties. Clearly, 
one needs a two-loop analysis when dealing with one-loop 
corrections of order 20%, in order to determine a s at the 
level of a few percent. Threshold resummation can also 
be performed. One should also take into account the 
one-loop electroweak corrections |6j. The QED ones are 
identical in form to the QCD ones, with the appropri- 
ate change in normalization, given by (0) . There will be 
a minor dependence on m t , which will be lessened with 
future Fermilab runs. The uncertainty on mt translates 
to an uncertainty on z. Since the j — cross section is 
near its peak for z ~ 1.2, minor variations in z will not 
appreciably affect the results. 

Of some concern are resolved photon contributions, 
where a gluon or quark within the photon can participate 
directly in the interaction. Suppression of these contri- 
butions is a major reason for working close to thresh- 
old. Since the parton distributions within the photon 
drop steeply with increasing momentum fraction, x, and 
since x must be large near threshold, such contributions 



are quite suppressed. Confirmation of this assertion may 
be inferred from the resolved contributions to b quark 
production near threshold presented in jl5j from which 
we conclude that only very poor knowledge (if any) of 
the photon structure will be required, as such contribu- 
tions will be a fraction of a percent of the cross section. 
One can further reduce those contributions by identify- 
ing outgoing jets collinear with one of the photon beams, 
which are a signature of resolved photon events. One 
can also require that the energy deposited in the detec- 
tors be equal to the total beam energy in order to account 
for missed jets of the type mentioned above. 

From the experimental side, we are assuming only 
that it events can be clearly identified. With experience 
gained from Fermilab, this seems reasonable, especially 
considering the cleaner initial and final states in the 77 
case. Another experimental issue is that of normaliza- 
tion. In order to avoid normalization uncertainties, aris- 
ing from luminosity uncertainties, we suggest the mea- 
surement of a ratio of cross sections denoted 



R 



77 _ ^(77 QQ + x) 



Q/P 



o-(77 



PP + X) 



P = W,l. (76) 



The ratio of it to W + W~ events is statistically quite 
powerful as over one million W + W~ events are expected 
at such a "W factory" jl(|. This highlights the comple- 
mentary nature of top quark and production at a 
77 collider. As well, electroweak corrections to W + W~ 
production have been studied |l6| . For the same reasons 
as for it production, the resolved photon contributions 
will be suppressed. If a bb pair is produced in conjunc- 
tion with the W + W~, this will constitute a background 
to it production. 

It is worth discussing the many advantages of deter- 
mining a s via 77 — > ti+X relative to some of the options 
currently being used. The calculation is perturbative and 
avoids nonperturbative contributions arising in a s deter- 
minations from mass splittings and tau decays. Other 
determinations, based on evolution of hadronic struc- 
ture functions, rely on the parton model and assumed 
knowledge of hadronic structure. No such assumptions 
are made here. Unlike the 3- to 2-jet ratio from e + e~ 
annihilation, wc avoid having to define the jet isolation 
criteria by measuring the total it cross section. Since we 
are at a large energy scale, not only does perturbation 
theory work well, but we automatically determine a s at 
(or above) the it threshold, without having to perform 
evolution or cross flavor thresholds. From a theoreti- 
cal viewpoint, the most comparably clean determination 
comes from the ratio of hadrons to lepton pairs produced 
in e + e~ annihilation at the Z pole. As mentioned ear- 
lier, the small QCD correction proves an insurmountable 
limiting factor in that case. 

At this stage, our enthusiasm is dampened somewhat 
however by the need for a two-loop calculation. This 
need is highlighted by the fact that there is an arbitrari- 
ness in the choice of renormalization scale, fi, which can 
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only be compensated by the inclusion of two-loop correc- 
tions. The variation of a s with In \i is order o? s though, 
so for a reasonable choice of [i (i.e. y/s, m t , ...) the 
two-loop scale dependent contribution should not be too 
large and should not change the value of a s radically. 
Nonetheless, as pointed out earlier, a two-loop calcula- 
tion will eventually be required. In light of that fact, we 
see the importance of having simple analytical results for 
the one-loop corrections as they will be incorporated in 
the two-loop result. Also, we do not suggest that one 
consider this determination of a s in isolation. Rather, it 
should be combined with all other precision determina- 
tions, including the low energy ones, in order to minimize 
the error and provide an excellent test of QCD at the 
same time. 



VIII. CONCLUSIONS 

The analytical results for the one-loop QCD/QED cor- 
rections to massive fermion production presented in |i| 
were extended by analytically integrating the single inte- 
gral (virtual+soft) part. The differential and integrated 
cross sections were series expanded to order /3 10 (includ- 
ing /3 11 ln/3 terms) and were shown to be of practical use 
as well as being informative. Accurate parametrizations 
of the total cross section, valid up to yfs/2m = 20, were 
presented. As an application, we showed how top quark 
production at a 77 collider capable of reaching y/s ~ 420 
GeV could be used to precisely determine a s , statistics 
permitting. Theoretical uncertainties were briefly dis- 
cussed as were advantages over other a s determinations. 
Those advantages make this method of determination 
quite appealing. 

The importance of performing the two- loop corrections 
was emphasized. As well, the best value of a s will still 
come from combining all determinations, including the 
low energy ones. The major strengths of this determi- 
nation are the largeness of the QCD corrections and the 
potential to reduce the theoretical systematic errors via 
inclusion of higher order contributions, calculable using 
well-established perturbation theory. The major weak- 
ness being the tedious nature of the required two-loop 
calculation. 
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APPENDIX A: 



Here we present the coefficients <Xj(j) appearing in the expression (^6[) for the analytically integrated single integral 
part, flp. They are 

01 = (4+14/3-6/3 2 )(l+/3)+j(7 - 9/3 - 19/3 2 + 3/3 3 ), a 2 = 8[2/3 2 (l + (3 2 ) - j(5 + 3/3 2 )], 

a 3 = -16 + 8/3 4 + j(6 - 14/3 2 ), a 4 = -8(1 - /3 4 + 2j)(l + /3 2 ), 

a 5 = -{[4(1 + 13 2 ) 2 - (3(1 - /3 2 )](1 - (3 2 ) + j[8 + (3/2 + 14/3 2 - 7/3 3 - 2/3 4 + 5/2/3 5 ]}, 

a 6 = (1 - (3 2 ) 2 - , ? /2(l - 14/3 2 + 5/3 4 ), a 7 = 6(1 - /3 2 ) 2 + j(25 - 3(3+ 10/3 2 + /3 3 - 9/3 4 ), 

a 8 = [(1 - (3 2 ) 2 + j(15/2 + 3/3 2 - 5/2/3 4 )]/2, a 9 = -(1 - /? 2 ) 2 + j/2(l - /3 2 )(17 - 5/3 2 ), 

oio = 4[(2 -(3 + 2/3 2 )(l + /3) + j(3 - 19(3 - 3(3 2 - 7/3 3 - 2/3 4 )/(3 + /3 2 )], 

an = [(1 + (3 2 ) 2 - 0/3(1 - /3 2 )](1 - /3 2 ) + j/6(12 + 29/3 + 39/3 2 - 10/3 3 - 9/3 4 - /3 5 ), 

012 = {[(1 + /? 2 )(9 + 9/3 + 9/3 2 + 2/3 4 ) - 3(1 - /3 2 ) - 2/3 3 ](l - (3)(3 

+j/2(180 - 219/3 + 165/3 2 - 219/3 3 + 95/3 4 - 73/3 5 + 35/3 6 - /3 7 + 5/3 8 )/(3 + /3 2 )}/(3 + /3 2 ), 

da - -[7(1 - /3 2 ) 2 + j/2(73 + 58/3 2 - 35/3 4 )]/2, a 14 = 8(1 + (3 2 - 7/2 j), 

Oi6 = -16[1 + (3 2 - 3j], oi8 = -4(1 + /3 2 )(1 - /3 4 + 2j), 

o X7 - 8[1 + 2(3 2 - j(19 + 7/3 2 )/(3 + /3 2 )], (Al) 

da = 4[(9 + 2(3 2 ){1 - (3 4 ) - j/2(33 - 27/? 2 - 29/3 4 - 9/3 6 )/(3 + /3 2 )]/(3 + (3 2 ), 

aw = [8(1 + /3 2 ) 2 - 3/3(1 - /3 2 )](1 - /3 2 ) + j(16 - 37/ 2/3 + 28/3 2 - 11/3 3 - 4/3 4 + 15 /2/3 5 ). 



APPENDIX B: 



Here we present the Cj entering in the parametrizations for /^(+,+), / (+> — ) i n t ne various energy regions 
whose form is given in Equations (74), dTq), respectively. For /^'' e (+, +), the coefficients are 



For /(!).«(+,+), 



For / (1) >' le (+,+), 



For /( x )' ie (- 



For /«>«(+,-), 



For / (1) ^ e (+,-), 



c 3 = -155, c 4 = -125.68, c 5 = 119.09, 
c 6 = -540.34, c 7 = 364.26 . 



c = 53.502, d = -137.56, c 2 = 102.57, c 3 = -29.359, 
c 4 = 3.3413, c 5 = .21711, c 6 = -.061446 . 



c = 71.912, ci = 47.622, c 2 = -.67576, 
c 3 = -2.1675 x 10~ 2 , c 4 = 1.1221 x 10 



-3 



c 5 = -667.218, c 6 = 2252, c 7 = -5395.85, 
c 8 = 8137.55, c 9 = -6654.48, c 10 = 2304.95 



c = 55.267, a = 57.115, c 2 = -7.3405, 
c 3 = 1.3777, c 4 = -.11197 . 



c = 207.66, a = 37.016, c 2 = -.15617, 
c 3 = -1.2899 x 10" 3 . 



(Bl) 



(B2) 
(B3) 

(B4) 

(B5) 

(B6) 
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FIG. 1. The functions (a) / (1) ; (b) / (0) , versus /?, in the 
low energy region, for the various helicity states. 
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FIG. 2. The functions (a) (b) / (()) , versus z, in the FIG. 3. The functions (a) (b) / (0) , versus z, in the 

intermediate energy region, for the various helicity states. high energy region, for the various helicity states. 




FIG. 4. The 77 —> tt + X cross section at LO and NLO, 
versus z, for the various helicity states. 
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TABLE I. The various for values of 1.2 < z < 20, and the corresponding single and double integral contributions. Here 
n.s. is the number of significant figures after the decimal point in and f£\ 



z 


n.s. 


/£> ('+>+) 




At) 

J unp 






/i 1)d (+°-) 


/«(+!-) 


1.2 


4 
4 


70.578894 
103.9951788 


-5.47998017 
-6.85659845 


65.0989 
97.1386 




33.4162848 
37.1626092 


-1.37661828 
-4.10336189 


32.0397 
33.0592 


2 


3 
3 


68.5516 
-7.8931 


-24.064 
14.728 


44.488 
6.835 




-76.4447 
144.9963 


38.792 
-62.856 


-37.653 
82.140 


3 


3 
3 


92.2075 
-99.5479 


-29.594 
96.259 


62.614 
-3.289 




-191.7554 
283.9629 


125.8526 
-155.447 


-65.903 
128.516 


4 


3 
3 


132.0495 
-153.7108 


-33.0381 
182.4102 


99.011 

28.699 




-285.7603 
417.8098 


215.4483 
-248.4864 


-70.312 
169.323 


o 


Q 
O 

3 


1 7fi 701 A 

-190.8526 


-ou.ouz 

263.100 


72.247 




Q«7 KKAf\ 
-OU I .OOIU 

544.2554 


9QQ Q09 

-336.704 


-u t .\jo£ 
207.551 


10 


3 
3 


395.3262 
-292.8618 


-55.3625 
584.182 


339.964 
291.320 




-688.1880 
1083.5142 


639.5445 
-694.907 


-48.6435 
388.607 


20 


1 
1 


749.8886 
-390.5080 


-79.437 
1007.530 


670.45 
617.02 




-1140.3966 
1890.2852 


1086.967 
-1166.404 


-53.43 
723.88 


TABLE II. The fractional errors 
1.05 < z < 1.4. 


on the various 


computed using 


the 


series expansions up to order f3 10 , 


for values of 


z 




P 


f.err(+,+) 


f. err(+, 


") 




f . 6IT un p 




1.05 
1.1 
1.2 
1.3 
1.4 




.3049 
.4166 
.5528 
.6390 
.6999 


2.1 x 10~ 7 
-6.8 x 10~ 6 
-2.0 x 10~ 4 
-1.4 x 10~ 3 
-5.6 x 10~ 3 


2.5 x 10 
2.3 x 10 
3.3 x 10 
1.3 x 10 
2.9 x 10 


-6 
-4 
-3 
-2 
-2 




4.3 x 10~ 7 
3.1 x 10~ 5 
6.9 x 10~ 4 

3.4 x 10~ 3 
8.9 x 10~ 3 


2.1 x 10~ 6 
6.6 x 10~ 5 
1.5 x 10~ 3 
7.3 x 10~ 3 
2.0 x 10~ 2 
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